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1 Introduction 



The explanation behind electroweak symmetry breaking (EWSB) is still largely unknown. As- 
suming the symmetry is spontaneously broken through the Higgs mechanism p^j still leaves 
room for speculation on the structure of the Higgs sector. Extending the standard model (SM) 
with softly broken supersymmetry at the TeV scale has several benefits: it protects the scalar 
masses from being sensitive to the large hierarchy of scales through quadratic divergences, it 
allows for gauge coupling unification, and it presents a candidate for the cold dark matter of 
the universe in terms of a stable lightest supersymmetric particle (LSP). 

In the minimal supersymmetric standard model (MSSM), the electroweak symmetry is 
broken by the presence of two complex Higgs doublets. The MSSM Higgs sector has a rich 
phenomenology which has been explored in exquisite detail [2]. Searches for MSSM Higgs 
bosons constitute a cornerstone of experimental particle physics at colliders. Driven in parallel 
by the negative results from direct searches, and indirect results from the measurement of 
e.g. electroweak precision observables, limits can be placed on their masses and couplings. The 
combined data implies that a SM-like Higgs boson should be just around the corner 0, and 
definitely within the reach of the LHC There are also theoretical arguments favouring a lightest 
Higgs mass rrih < 130 GeV in the MSSM, where < mz hold at tree-level. It has even been 
argued [1], that the MSSM gives a better combined fit to the available data than does the SM, 
adding further support to the idea of weak scale supersymmetry. 

Despite the great interest in the minimal model, we find it advisable to keep an open mind 
also with regards to further enlargements of the Higgs sector. After all, there exists no direct 
evidence for any particular theory beyond the standard model. It must therefore be ensured 
that no experimental opportunities are lost due to ignorance of parameter regions which are 
excluded in the MSSM. Corresponding parameter regions may be allowed in other models, 
as we shall demonstrate explicitly in this work for the simplest extension of the MSSM, the 
next-to-minimal supersymmetric standard model (NMSSM). For a general introduction to this 
model, we recommend the recent reviews O [6] . Details necessary for this paper, in particular 
on the extended Higgs sector, will be described below. 

Let us briefly recall the main theoretical motivation for the NMSSM. The MSSM superpo- 
tential W^*^^ contains, in addition to the Yukawa terms of mass dimension d = 3, the bilinear 
coupling H^MSSM — l^Hu-Hd of the two Higgs doublet superfieldsj^ The parameter /i has positive 
mass dimension, but since it is introduced directly in the supersymmetric theory it is not associ- 
ated with any obvious scale (e.g. the SUSY breaking scale). This leaves only /x = or Mqut as 
natural choices, both of which are phenomenologically impossible. An acceptable phenomenol- 
ogy requires /i close to the scale of electroweak symmetry breaking: ~ /xew < ^susy, 
where Msusy is the scale of supersymmetry breaking. This is known as the '/i problem' of the 
MSSM. The primary motivation for the NMSSM is to solve this problem. Replacing Wj^ggjyf 
with a coupling of the Higgs doublets to a new singlet field 5", the /x parameter is generated 
dynamically from soft terms as the vacuum expectation value of 5*. The NMSSM can also help 
to reduce the required amount of fine-tuning compared to the MSSM 

^We use the notation where $ denotes a chiral superfield and <& its scalar component. The dot gives the 
S'C/(2)-invariant scalar product such that ■ Hd = H:^H^ — H^H^^. 
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The NMSSM does not only have favourable properties from a theoretical perspective, but 
it also has an interesting Higgs phenomenology which extends that of the MSSM. The singlet 
components mix with the Higgs scalars (and neutralinos) , modifying the tree-level mass rela- 
tions, and introducing new modes of Higgs production and decay which can be relevant for the 
LHC. The present work starts from the MSSM limit of the NMSSM, which will be discussed 
further in the next section. More specifically, we consider low-energy MSSM benchmark scenar- 
ios as our starting point. These are easy to work with due to the relatively low number of free 
parameters; in the Higgs sector typically the CP-odd Higgs mass (or equivalently mH±) and 
tan /3. The simplicity inherent in the benchmarking approach is appreciated by experimental 
collaborations, sometimes to the degree that one could fear loss of information in the process 
of reporting results. One of the aims of this paper is therefore to use present experimental 
limits to determine how our current understanding of well-known scenarios in the MSSM is 
modified when going to the NMSSM. To this end, we use input from LEP, the Tevatron, and 
flavour physics experiments performed at lower energies to constrain the parameter space. Of 
special interest to us is the possibility to have comparably light Higgs bosons, while still being 
compatible with all known limits. This question is investigated in more detail, focusing in par- 
ticular on the cases where a CP-odd Higgs and/or the charged Higgs boson is very light. For 
scenarios that fulfil the experimental constraints, we calculate and discuss Higgs decay modes 
and production cross sections for LHC running at 14 TeV. 

This work does not concern the idea of a hidden NMSSM Higgs sector, or how to establish a 
'no- loose' theorem under these circumstances [HI E] . Quite to the contrary, we focus on scenarios 
with one or more light Higgs bosons, maximizing the accessible phenomenology while avoiding 
exclusion by the present data. Work along similar lines is presented in [10] , with the motivation 
of explaining the slight excess of Higgs-like events at LEP observed around rrih — 90 GeV as 
Hi production followed by the dominant (but unobserved) decay Hi — )■ AiAi [TH [T2]. The 
phenomenology of these so-called ideal Higgs scenarios — although obtained through a different 
strategy — should be closely related to some of the scenarios we investigate here. In addition to 
our different starting point, we extend the earlier results by including new limits from Higgs 
searches at the Tevatron in our analysis, as well as the charged Higgs boson searches at LEP 
and the Tevatron. We also present combined results from collider and flavour constraints. Our 
independent implementation of the latter serves as an important cross-check and an update on 
the Higgs constraints in the NMSSM. Finally, we show explicitly the effects of these constraints 
on the quantities of phenomenological interest at colliders. Unlike the NMSSM benchmark 
scenarios presented in p3] — where the charged Higgs boson is always taken to be heavy — we 
also consider the case when it is light; an interesting possibility discussed in [T?]. 

This paper is organized as follows: in Section [2] we introduce the NMSSM Higgs sector and 
describe briefly the computation of the Higgs masses and mixing at tree-level. This is followed 
in Section |3] by a more detailed description of our analysis strategy. Section |4] first presents the 
experimental input, which is then applied to give up-to-date exclusion limits on the NMSSM 
Higgs boson masses. The LHC phenomenology for the Higgs sector of models allowed by the 
combined constraints is discussed in Section |5| Finally, Section [6] contains a summary and the 
conclusions. 
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2 The Higgs sector of the NMSSM 



To simplify the analysis, a common restriction of the most general supersymmetric model with 
a gauge singlet is to consider a scale invariant superpotential. This model is sometimes referred 
to as the Z3-symmetric NMSSM. With this assumption, the superpotential takes the form 

WnMSSM = W^MSSM + XSHu ■ Hd + g'^'S'^, (1) 

with all parameters dimensionless. W^mssm contains the Yukawa terms of the MSSM which are 
not modified. The supersymmetric parts of the Higgs potential are derived by evaluating the 
usual F- and D-terms. In addition to the superpotential, a complete phenomenological model 
requires specification of the soft SUSY-breaking potential to have Vniggs = Vp + Vd + Vjoft- The 
terms in Vson containing only the Higgs fields are given by 

Vsoft = m]jJ^Hu\^ + m]j^\Hd\'' + ml\S\^ 

+ {XAxSH^ ■ Ha + ^kA^S' + h.c.) . (2) 

The remaining terms in Vgoft keep their MSSM expressions. In addition to the unstable dou- 
blet minimum required for EWSB, it is assumed that m| < Af./9 (which may be radiatively 
generated) to make S = unstable. Expanding the singlet field "S* = f s + {4>s + around 
its vacuum expectation value Vs = {S) gives rise to an effective parameter /Xgfr = Xvs with 
mass dimension. Since Vg originates from SUSY breaking, a phenomenologically viable value 
|/ieff| ^ ^susY is natural. Having noted this important dynamical mechanism to generate the 
H parameter, we drop the 'eff ' subscript in the following. 

Following EWSB, the physical spectrum contains seven Higgs scalars, five of which are 
neutral, and a single charged pair. If CV is conservecj^ (as assumed here) the neutral mass 
eigenstates are arranged into three CP-even Higgs bosons {Hi, H2, H^), ordered in mass such 
that rriH-^ < 112^2 < f^Hs, and two CP-odd states {Ai, A2) with m^^ < m^ij. The charged Higgs 
boson — which makes a proper mass eigenstate also in the presence of CV violation — is denoted 
by H^. Using the minimization conditions for the potential, the scalar mass parameters of 
Equation ^ can be eliminated in favour of the vacuum expectation values of the doublets. We 
define as usual 

vl + vj = v^ c^iUAGeVY, tan/3 = 

Due to the additional interactions with the NMSSM singlet that have been introduced, a full 
specification of the Higgs sector now requires six parameters at tree-level: 

X K Ax Ak tan/3 Vg. 

Conventions can be chosen such that A and tan /3 are positive, while the remaining parameters 
can have either sign. 

^Unlike the MSSM, automatic CV conservation at tree-level is not guaranteed, but constitutes an additional 
assumption. In the CV conserving case it is generically possible to use real values for all parameters appearing 
in the Higgs potential. 
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The presentation of the Higgs mass matrices at tree-level is simplified by defining an effective 
doublet mass rriA as 

m\ = — -(Ax + kVs), (3) 
S/3C/3 V / 

with Cjs = cos (3, and Sjs = sin (3. Note that in general does not correspond to the mass 
of any physical Higgs boson in the NMSSM. Eliminating the massless degree of freedom eaten 
by the Z, the remaining 2x2 mass matrix for the pseudoscalars can be written in the basis 
{A, as), where A = Im(cos f3Hu + sin is the doublet component. In this basis, it is given 

by 

Ml = 

V v'^c s 

— (m^c^S/3 - 3XkvI) ^ {m\ci3Si3 + 3A/tt;^) - 3kA,,Vs , 

The pseudoscalar masses are diagonalized by a simple rotation with an angle 6 a such that 

Ai = A cos 6a + o-s sin 6a / , x 

A2 = —As'm6A + crsCos6A, 

where an explicit relation for the mixing angle is 

cos6'^ = — ' (5) 



-^l>,12+«--^P,ll)" 

To describe the mixing in the CP-even Higgs sector requires the introduction of a full 3x3 
unitary matrix. In a basis 5"^'^'^'^ = (Re(ifu), Re{Hd), Re(S')) where the mass eigenstates Hi are 
given by Hi = SijSj"^^^ we use the parametrisation 

(C12C13 S12C13 Si3 

— ■512C23 — C12S13S23 C12C23 — S12S13S23 C13S23 I , (6) 
— C12C23S13 + S12S23 — C12S23 — S12S13C23 Ci3C23^ 

with Cij = cos 6ij and Sij = sin 6ij^ For brevity, we do not give the entries of the CP-even 
mass matrix Ais here. They can be found for example in [U]. The tree-level value of rriHi is 
no longer limited from above by < m|cos^2/3 (as is the case in the MSSM). Instead the 
modified limit 1151 



m 



2^2 
w 



j,^ < cos^ 2/3 + sin' 2(3 (7) 



applies. Finally, the tree-level mass of the physical charged Higgs boson is given by 

2 2,2 \2 2 

m^± = + m^y - A f , 



■^This corresponds to a mixing in the MSSM limit which takes the form 6*12 — > —a, cosOi^ — > 1, and 

cos 6*23 1- 



5 



where m\ is given by Equation (|3|. For a fixed doublet mass, the term — A^w^ reduces the 
charged Higgs mass compared to its MSSM value. 

The whole discussion so far concerned the Higgs masses at tree-level. Like in the MSSM, 
the tree-level relations may be subject to sizeable higher order corrections [121 [I7|. The Higgs 
mass corrections taken into account in our numerical analysis are discussed further in the 
next section. As a final point, let us discuss the MSSM limit. This is obtained by taking 
simultaneously A — i- and k — )■ 0, keeping the ratio k/X constant, thus decoupling the singlet 
from interacting with the Higgs doublets. A finite > in this limit (as required to get 
massive charginos) is achieved by taking — )■ oo. The remaining dimensionful parameters are 
held constant. This illustrates the possibility to consider a class of models with a continuous 
transition from the MSSM to the NMSSM. Finally, we note the interesting fact that even if 
all influence on Higgs physics vanishes the LSP may be singlino-like, giving a possibly long- 
lived (or even charged) NLSP which could lead to modified collider phenomenology even in the 
decoupling limit [IB] ■ 

3 Analysis Strategy 

The masses (and mixing matrices) of the Higgs bosons receive numerically important corrections 
at higher orders. These loop corrections introduce a dependence in the Higgs sector on the 
soft SUSY-breaking parameters of other sectors, most notably from stops (for high tan /3 also 
sbottoms), which have sizeable couplings to the Higgs doublets. Both the dominant corrections 
proportional to the Yukawa couplings yf, yl and the mixed electroweak/ Yukawa terms of order 
d'^Uti d^Vb similar between the MSSM and the NMSSM. For the numerical analysis of the 
NMSSM Higgs sector we use the code NMSSMTools 2.2.0 [19J, which includes the dominant 
one-loop and leading logarithmic two-loop corrections to the Higgs masses. 

To study how constraints on the Higgs sector are modified when going from the MSSM to 
the NMSSM, we adopt a strategy to fix the higher order corrections from sparticles in order to 
focus on 'genuine' NMSSM effects which appear already at tree-level. To achieve this, we make 
use of MSSM benchmark scenarios for Higgs searches which were originally presented for LEP 
|20j . and later extended to hadron colliders [21]. The benchmark scenarios provide a starting 
point in terms of the soft SUSY-breaking parameters at the electroweak scale. Recall that the 
parameters specified in these scenarios are the universal scalar mass Msusy, the gaugino masses 
M2 and M3, with Mi typically given by the GUT relation Mi = ^^^^2, the trilinear couplings 
for the third generation At = = Ar, and the value for /i. Once everything is fixed in this 
way, only rriA (or mH±) and tan/3 remain as free parameters. 

Since we are interested in Higgs mass limits, we look at benchmark scenarios covering the 
two extremes of radiative corrections to the lightest Higgs mass: one scenario with maximal 
mixing in the stop sector and one with no mixing. The full parameter sets for the two scenarios 
we consider are given in Table Q Scenario (A) is the constrained maximal mixing scenario 
and scenario (B) is called no mixing. Mixing refers to the size of the off-diagonal entries in 
the sfermion mass matrices, which for the stops are given by rritXt where Xt = At — /icot /3. 
When the mixing is maximal, it also maximizes the one- loop corrections to the lightest Higgs 
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Parameter 


Scenario (A) 


Scenario (B) 


^SUSY 


1 TeV 


1 TeV 




-\/6MsusY 







200 GeV 


200 GeV 


M2 


200 GeV 


200 GeV 


M3 


800 GeV 


800 GeV 



Table 1: Values for the soft SUSY-breaking parameters and /i in the MSSM benchmark scenarios 
(A) constrained maximal mixing and (B) no mixing of [21] • In the on-shell renormalization 
scheme for rriA, iXjr'^l = 2Msusy for maximal mixing. 

mass. Including also two-loop corrections, the value of Xt for which this is achieved depends 
on the renormalization scheme; see [22] for a discussion and comparison between MS and OS 
(on-shell) results. It has been shown that certain two-loop terms lift the degeneracy of the 
ruh maximum under Xt — )■ —Xt. This is of some significance to us here, since we consider the 
constrained maximal mixing scenario which differs from maximal mixing only in the sign of 
Xt. Even though this gives a lower maximum value for tuhi (by up to 5 GeV), this scenario 
is in better agreement with data from 6 — )■ S7 transitions]^ A negative Xt at the weak scale is 
also strongly favoured in constrained MSSM scenarios with parameters evolved from the scale 
of grand unification. For the scenario with no mixing, Xt = Xb = X^ = is imposed. To 
calculate the MSSM Higgs spectrum we use the code FeynHiggs 2.7.4 [23] . 

Our strategy for the NMSSM case is to extend the MSSM benchmarking approach in the 
following manner: we keep the soft parameters fixed at the values given by the MSSM scenario. 
Also the effective value of /i = Xvg is kept at its MSSM value. As before, this leaves the 
tree-level Higgs sector free. Since the CP-odd Higgs mass parameter is not physical in the 
NMSSM, we would like to use mj^± as input. This is achieved by using Equations ^ and (|8| to 
determine Ax at tree- level, which is then corrected at higher orders (using NMSSMTools) through 
iteration. In summary, the NMSSM scenarios have three additional parameters compared to 
the MSSM: 

{mH± , tan (3) — )■ {mH± , tan /3, A, k, A^) . (9) 

To present constraints on the NMSSM Higgs sector in a comprehensible way, we will perform 
scans over the additional parameters. Their domain of interest is limited from the properties of 
the theory. Following [21], we require that A, and the Yukawa couplings (most importantly 
yt, or for large tan/3 also Hb) remain perturbative to the GUT scale. This leads to the upper 
limit 

VA2 + ^2 < 0.7. (10) 

This condition is illustrated in Figure [l| which shows for scenario (A) a sample of weak- 
scale values for (A, k) compatible with perturbativity. The same analysis gives a lower limit on 
tan/3 > 1.4 using the top mass rrit = 173.3 GeV [25J. We employ this value for the top quark 
pole mass throughout this work. 

■*The main sparticle contribution to BR(_B — >■ Xgj) is proportional to sgn(/iAf). Since the other major 
contribution (from H^) is always positive, the negative sign for At is required for their successful cancellation. 
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Figure 1: Distribution of allowed points in the space of NMSSM couplings for scenario (A). 
Dark (black) points are compatible with perturbative unification, while the brighter (green) 
points in addition fulfil collider constraints on Higgs masses. 



The range for can be restricted by demanding the symmetry breaking minimum of the 
potential to be stable. From this can be derived [23] an approximate condition valid for high 
tan /3 and high 

-'f<A.<0. (11) 

Since this condition is only approximate, and since positive scalar mass squares are necessary 
for stability, we do not employ scanning over yl^ directly. Instead we use the lightest CP-odd 
Higgs mass rriAi as input parameter and determine through an iterative procedure. The 
stability of the EWSB minimum is checked numerically. Since the charged Higgs mass does 
not depend on A^ to the computed order, we can use both these masses as independent input 
parameters in a pseudo on-shell approach. We find that this greatly improves scanning efficiency 
for low Higgs masses, with about 60% of the tried points giving a physically viable spectra for 
the scanning ranges described below. The distribution of points in the {ruAi 5 ''t^h^ ) plane 
obtained for scenario (A) using this method is shown in Figure |2} 

It is important to exercise caution when interpreting results from scans over multi-dimensional 
parameter space, especially when the sampled parameters are related non-linearly to the funda- 
mental parameters of the theory as is the case here. In particular one should not interpret the 
posterior density of points in a particular region as a probability in the frequentist sense. What 
can be shown by the procedure employed here is the existence of valid models or — trusting the 
sampling is sufficiently covering — their non-existence in absence of valid points for a certain 
parameter region. With this objective in mind, we perform random scans over the parameter 
space of each benchmark scenario with uniform sampling in the ranges mH± & [80,500] GeV, 
tan/3 G [1.2,50], A G [0,0.6], k G [-0.7,0.7] (requiring Equation ^ to hold), and finally 
G [4, 500] GeV with the upper limit < imposed by the tree-level relation given in 
Equation ([s]). We limit ourselves to the case m^i > 2m^ (where hadronic decay modes other 
than bb are not relevant). Constraints on an even lighter Ai are presented in |26j. For each 
scenario, we generate (9(10^) physical parameter points. 
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Figure 2: Masses of the lightest CP-odd Higgs boson Ai and charged Higgs boson in the 
NMSSM for scenario (A). All accepted points in the displayed mass range are shown in black, 
while points that are in addition compatible with collider constraints (see Section |4]) are shown 
in green. 

4 Experimental Constraints 

To investigate the viability of the considered scenarios, we take into account existing constraints 
from several different sources as detailed below. As a general constraint, the direct limits on 
sparticle (not Higgs) masses compiled by the PDG |27] are checked. In the benchmark scenar- 
ios considered here, they are always fulfilled. Furthermore, to provide a possible explanation 
for the dark matter density of the universe, it should be required that the lightest supersym- 
metric particle (LSP) is neutral. Here this is the lightest neutralino which has a typical mass 
m^o ~ 90 GeV. Since the density Q^^h"^ of cold dark matter (CDM) in the universe is very 
precisely determined by WMAP [28] in the so-called concordance model, strong constraints 
can be derived on the SUSY parameter space under standard thermal relic assumptions for 
the nature of the DM. This constraint can be used to reduce the effective dimensionality of 
parameter space by one unit. However, we do not apply the CDM constraints here, since we 
are interested in the Higgs sector, and these constraints could equally well be fulfilled by mak- 
ing a change to the underlying MSSM parameters (similarly to |29]). Moreover, it is always 
possible to accommodate the WMAP constraints by changing the cosmological assumptions in 
the early Universe (yet compatible with the cosmological observations) [30]. A detailed study 
of this topic goes beyond the scope of this paper. 

4.1 Collider constraints 
4.1.1 LEP 

The four LEP experiments have placed limits on the masses and couplings of neutral MSSM 
Higgs bosons in a large number of channels |31], most of which are also relevant in the NMSSM 
case. Primarily two modes of neutral Higgs production were probed at LEP: the Higgs-strahlung 
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process e"'"e^ — t- ZHi, by which a single Higgs boson is produced for each event, and the 
associated production e"'"e~ — )■ HiAj. The Higgs decay modes considered in the LEP analyses 
involve most final state combinations of rif, bottom quarks and n — nf, tau leptons for n,nf, < 6. 
There also exists an important limit on the lightest CP-even Higgs mass rrih in the MSSM 
from a decay mode independent search where the mass spectrum of the recoiling Z boson was 
analysed. We take into account the exclusion by all these channels (and many others, in total 
541 channels for NMSSM case) consistently at 95% CL using the code HiggsBounds 2.0.0 
|32j which contains a collection of data on Higgs searches]^ One particular channel — not part 
of the original LEP analysis and therefore not yet present in HiggsBounds — is Higgs-strahlung 
followed by the decay Hi — )■ AiAi — )■ t~^t~t~^t~ for a very light Ai {rriAi < 2mf,). This 
mode received recent attention as a potential probe for light pseudoscalar Higgs bosons in the 
NMSSM [12], which spurred a reanalysis of old data and the publication of first experimental 
results in the interesting mass range by ALEPH [33j. These new results are included in our 
analysis. 

In addition to the neutral Higgs searches, the LEP experiments also published negative 
searches for charged Higgs bosons in the pair production channel e~^e~ — )■ H~^H~ , followed by 
the decays — t- r^Ur, — )■ cs, or — > W^A in all the different combinations [31]. The 
latter mode is not relevant in the MSSM, but could become important in the NMSSM (or a 
general two-Higgs-doublet model) with a sufficiently light Ai. The resulting direct limits on 
mi{± depend slightly on the precise branching ratios, but in general the combined lower bound 
is around mH± > 80-90 GeV. We use mH± > fnw as the lower limit for our NMSSM scans in 
this study. The precise charged Higgs boson limits are taken into account using HiggsBounds. 

4.1.2 The Tevatron 

The two Tevatron experiments CDF and D0 are actively searching for Higgs bosons in pp 
collisions at 1.96 TeV, and their combined results currently exclude a SM-like Higgs at 2 a 
in the mass interval 158 GeV < < 175 GeV [3S]. This range is already beyond the 
theoretically accessible region in the models we consider, but for supersymmetric Higgs bosons 
with enhanced couplings to fermions, new modes that are currently not excluding for the SM 
Higgs become interesting. In particular the search for A production, followed by the decay 
A —7- r+r^, provides strong exclusion for high tan /3 due to the enhancement of the A coupling 
to down-type fermions. 

In addition to the neutral Higgs results, mass limits exist for a charged Higgs boson decaying 
in the — r^z/^ and — )■ cs modes from CDF [36] and D0 [37]. In the case of a 
light charged Higgs boson {mH± < rrit) the limits are reported on BR(t — )■ bH^), assuming 
either of the two modes saturates the width. In the NMSSM with a light Ai — where 
— )■ iy^y4i can be relevant — the experimental results can be conservatively interpreted as 
upper limits on the rescaled quantities ^h^cs = BR(t — )■ bH"^) x BR^H"^ — )■ cs) and ^h^tu = 
BR(t — )■ bH^) X BR(if^ — )■ t^Ut)- An additional constraint can be derived from the search for 
— )■ VT^Ai, Ai — ^ r+r^. This channel was recently analysed for the first time [38] . 

^The LEP limits are also available in NMSSMTools, and we have verified that applying these give results 
consistent with those from HiggsBounds within expected uncertainties. 
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Results from both neutral and charged Higgs bosons searches in the main channels at the 
Tevatron are incorporated in the HiggsBounds code. We use this to evaluate constraints both 
for the MSSM and the NMSSM. The charged Higgs constraints are cross-checked with an 
independent implementation where we also include the — )■ ly^Ai mode. 

4.1.3 Combined results 

We combine MSSM and NMSSM constraints from LEP and the Tevatron in each figure. This 
is demonstrated in Figure |3| which shows the lightest CP-even Higgs mass versus tan /3. The 
solid (coloured) regions are results obtained for the underlying MSSM scenario: the red (blue) 
areas are excluded by LEP (Tevatron) at 95% CL, while the green regions are allowed. White 
areas are theoretically inaccessible in the MSSM scenario. The scattered points (in black) are 
superimposed on the MSSM regions. They are the result of the NMSSM scans, with each point 
representing a choice of parameters compatible with all the collider constraints. 

Considering first scenario (A) (left two plots), a number of interesting features can be seen. 
There is very good agreement for a wide range of tan /3 between the upper limits on ntH^ in the 
MSSM (shown in the top row) and the NMSSM (lower). This demonstrates that the dominant 
radiative corrections are taken into account consistently, and that they indeed are similar. The 
upper limit on is higher in the NMSSM in a narrow band at low tan (3 ~ 2.5, where 
the tree-level niH^ is maximized by additional NMSSM contributions [39]. Focusing now on 
the excluded regions, it is clear from Figure |3] that the MSSM region of tan {3 < 10, where 
essentially the SM limit m/j > 114 GeV applies, is reopened in the NMSSM. Additionally, the 
almost tan /3-independent limit of nih > 90 GeV is not respected, but valid NMSSM points are 
observed down to mn^ — 60 GeV. That such light Higgs bosons avoid the collider constraints 
is due to a large singlet component of ifi, making its coupling to vector bosons smaller. In a 
situation where the H2 has a reduced coupling with respect to the SM C{H2VVY — 1, while 
the lightest Higgs boson has C^HiVV)"^ ^ 0, the lightest state Hi is not produced through 
Higgs-strahlung and the LEP limits do not apply. (For a precise definition of the reduced 
couplings we refer to Table [s] below.) 

The situation in scenario (B) (right plot) appears somewhat different, mainly because of the 
modified MSSM conditions. Close to the maximum value m/j ~ 115 GeV, only a narrow strip 
where tan /3 > 15 is allowed in the MSSM. On the other hand, there is still an allowed region for 
somewhat lower m^^^, ending at the decay mode independent LEP limit mh > 90 GeV, which 
remains similar to before. Even if the valid NMSSM points tend to cluster at large values for 
TTiHi, their existence in all the MSSM-excluded regions is evident. The region 110 GeV < < 
115 GeV is allowed because the decay Hi AiAi ^ 46 is open, and the LEP constraints for 
this mode are weaker than for the Hi — )■ bb decay which is always dominant in the MSSM for 
this mass range. 

For the lightest CP-odd Higgs boson Ai, exclusion plots of the same type as for Hi is shown 
in Figure |4j In scenario (A), a similar result as for Hi with respect to the LEP constraints 
can be observed: the MSSM limit rriA ^ 90 GeV is not respected by niAi in the NMSSM. The 
90 GeV limit applies only to the effective doublet mass (and then only when m^^ > 2mf,). 
We find points with a large doublet component also below the bb threshold. This supports the 
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Figure 3: Combined collider constraints in the (m//^,tan/3) plane for scenario (A) (left), and 
scenario (B) (right). The coloured regions represent exclusion in the MSSM at 95% CL from 
LEP (red), the Tevatron (blue), and allowed parameter space (green). The black points repre- 
sent NMSSM models which are compatible with all the collider constraints at 95% CL. 



results of [TU] , and shows that scenarios like this are viable with the updated collider constraints 
we have included. The most important constraint in this mass range is the new result from 
ALEPH on Hi — )■ AiAi — > 4r. It relates a light Ai to specific requirements on the mass of 
Hi. Either niHi > 107 GeV and the constraint does not apply, or H2 must be SM-like so 
that the coupling C{HiVV) of Hi to vector bosons is suppressed. In the intermediate range 
2mfe < niAi < 90 GeV, the doublet mass tjia has to be above 90 GeV. These values for tjiai 
are allowed in the NMSSM since the reduced coupling C{HiAiZ) can be lower than in the 
MSSM, which decreases the rate for associated HiAi production. We find that C{HiAiZy < 
0.5 is sufficient to find allowed rriAi ^^11 range (see also Figure 12 below). Having 

mAi <^ 90 GeV of course opens interesting phenomenological possibilities, as will be discussed 
further in Section [5] Another feature which is clearly visible in Figure |4] is the Tevatron MSSM 
exclusion at high tan /3 and intermediate rriA- This exclusion does not remain valid in the 
NMSSM when the relevant couplings C{Aibb), C{Ait^t^) oc tan (3 cos 9a are suppressed by 
singlet mixing (| cos^^il < 1). 

In scenario (B), the allowed NMSSM points divide into distinct classes. For tuai < '^^rrih, 
the same possibility as discussed for scenario (A) exists to have an Ai with a large doublet 
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Figure 4: Combined collider constraints in the (myij,tan/3) plane for scenario (A) (left), and 
scenario (B) (right). The colour coding is the same as for Figure IS] 



component. The same conditions on the mass of Hi as discussed for scenario (A) apply also 
in this case. In the region 2mf, < itia^ < 57 GeV we find points which are associated with 
a dominant decay Hi — )■ AiAi. As discussed above, a substantial branching fraction for this 
mode circumvents limits on rriHi, which leads to the large number of allowed points in this 
mass range. The upper limit is set by the kinematics from the upper limit on rriHi in scenario 
(B). In addition to these two classes, there is a sparse population of scattered points over 
most of the plane, preferentially in the MSSM-allowed areas. Where necessary, these avoid 
the collider limits by a small coupling C{HiAiZ), as also discussed for scenario (A). We note 
that tan /3 < 10 is disfavoured for > 57 GeV also in the NMSSM, since in this region it 
is difficult to find allowed values for rriHi- The only exception is for tan/3 ~ 2.5 where the 
tree-level value for rriHi is maximal. 

Figure [5] shows the combined constraints in the plane of the charged Higgs mass and tan/3. 
In the MSSM, these are quite similar to the results for the CP-odd scalar. This is a consequence 
of the tree- level mass relation = m^ + m^, which typically receives only small corrections 
at higher orders. As a result, an indirect limit > 120 GeV can be derived from the previous 
result rriA > 90 GeV. Turning now to the NMSSM results, it can be seen that this limit does 
not apply. There are two reasons why it is not effective: first — as we have already seen — even 
a doublet-like Ai can be very light (m^^ < 2m{,, which can correspond to an effective doublet 
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Figure 5: Combined collider constraints in the (mj:/±,tan/3) plane for scenario (A) (left), and 
scenario (B) (right). The colour coding is the same as for Figure [s] Constraints from direct 
searches for , as well as indirect constraints from neutral Higgs searches, are included. 



mass niA < 90 GeV). Second, and equally important, the tree-level mass relation is modified, 
cf. Equation ([s]). For this reason, we find the largest number of points with light if^ for large 
values of the NMSSM coupling A. Paradoxically, the MSSM limit (A — 0) is thus most clearly 
visible in the charged Higgs constraints — even though no new charged degrees of freedom have 
been introduced. In addition to the indirect constraints from the neutral Higgs searches, which 
are most excluding, mH± is limited by the direct search results. Here the decay — )■ W^Ai 
can degrade the standard search channels {H^ — )■ r^z/^-, — )■ cs), which would also make it 
easier to accommodate a light charged Higgs. However, we find that with the present (limited) 
reach of the Tevatron, it is not necessary to invoke this mode to avoid exclusion, save for very 
high (or very low) tan p. We expect this to change quite dramatically when results from the 
LHC become available. As a final remark on Figure [5} we would like to point out that the 
low tan P region is not excluded in the NMSSM. It could therefore be potentially interesting 
to search for two jet decays of H^. The situation for the charged Higgs boson in scenario (B) 
is very similar to scenario (A). Also in this case there exist points with 'mH± ~ 90 GeV which 
are still allowed. We note that a low tan/3 is only possible when the Ai is comparably light, 
cf. Figure |4} which means the decay — )■ W^Ai is likely to be dominant for tan /3 < 10 in 
this scenario. 



4.2 Flavour physics 

In addition to the collider constraints, we consider limits resulting from several flavour physics 
processes on the two benchmark scenarios (A) and (B). The same observables have been applied 
previously to constrain the Higgs sector in the MSSM in the context of the CMSSM and 
NUHM frameworks [IDl SD 112], and a subset of them also in the NMSSM [13]. For the 
theoretical evaluation of the model predictions, we use the code Superlso v3.0 [SI US] for 
the calculations both in the MSSM and the NMSSM. Superlso provides automatic interfaces 
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Observable 



Experimental value 



95% CL range 



BR(5 ^ X,7) 



(3.55 ±0.24 ±0.09) x 10"^ 



[48J 
[49J 
[48J 
[50J 
[51J 
[48J 



[2.16,4.93] X 10"^ 
< 4.7 X 10-^ 
[0.56,2.70] 
[0.151,0.681] 
[0.990,1.018] 



BR{B DTUr)/BR{B Deu^) 



< 4.3 X 10-8 
1.63 ±0.54 



Rm{K -> /iz/) 



0.416 ±0.117 ±0.052 
1.004 ±0.007 
(5.38 ±0.32) X 10-2 



[4.7,6.1] X 10-2 



Table 2: Experimental values for flavour physics observables and the corresponding allowed 
ranges at 95% confidence level. 

with several MSSM and NMSSM spectrum generators. Here we used SOFTSUSY 3.1.6 and 
SuSpect 2.41 [47] for the MSSM and NMSSMTools 2 . 2 . [I9] for the NMSSM. 

To determine the allowed intervals for the model to be compatible with the data at 95% CL, 
in the presence of both experimental and theoretical uncertainties, a Monte Carlo approach is 
used [52]. The experimental values used as input for our analysis and the resulting allowed 
ranges are given in Table [2j Since the main objective is to compare the two models, no attempt 
is made at a sophisticated statistical combination of the different constraints or a global fit. 
Instead we simply give the exclusion at 95% CL by each observable separately, meaning that 
any particular point may be excluded by more than one measurement. This approach should 
be sufficient to give a general indication of what the viable regions for the Higgs masses in the 
NMSSM are. 

The NMSSM specific additions to BR(B — )• Xs7) required to take into account the additional 
contributions from the Higgs bosons and the extra neutralino are obtained according to a 
generalisation of [531 El]- The NMSSM effects in this decay appear only at the two- loop level 
and therefore do not lead to substantial departure from the MSSM results. 

The implementation of BR{Bs — )■ is based on [Ml EHl ES] and is in the approximation 

of large tan/?. Contrary to BR{B — t- X57), here one can expect some differences with respect 
to the MSSM. Indeed, when Ai is light (which can be possible for any m/^±), it generates 
sizeable contributions to the branching ratio through penguin diagrams, especially at large 
tan/?. On the other hand, even for large itiai coupling of Ai with down type quarks can 
take both small and large values [56], and if this coupling is large the Ai contribution can 
become competitive or even dominant. For the numerical evaluations we use a value for the Bg 
decay constant of /b, = 238.8 ± 9.5 MeV [57]. 

In the leptonic (and semi-leptonic) decays of pseudoscalar mesons, the dominating con- 
tribution beyond the SM is through charged Higgs exchange diagrams. The calculation of 
these observables have been simply extended from the usual MSSM results to the NMSSM. 
Since this sector is not directly modified in the NMSSM, the impact on fiavour physics is 
only indirect. Small differences can arise due to the fact that the charged Higgs mass in the 
NMSSM can be lower than in the MSSM for a similar neutral Higgs spectrum. Here we con- 
sider the branching ratios of S — )• rz/^, B DrUr, Dg — )■ tu^ and K — )■ /iz/, and we use 
[Kfel = (3.92 ± 0.09 ± 0.45) x 10-^ [58], /b = 192.8 ± 9.9 MeV [57], /d, = 248 ± 2.5 MeV [59], 
and JkI fit = 1-189 ± 0.007 [20] for the numerical evaluations. 
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Figure 6: Flavour constraints on the (m^fi , tan /3) plane in the MSSM (upper) and the NMSSM 
(lower) for scenario (A) (left) and scenario (B) (right). The constraints are applied in the order 
indicated by the legend with the points satisfying all constraints in green in the foreground. 



All the formulae and the detail of the calculations of the above branching ratios can be 
found in |44j . 

In addition to the aforementioned observables, we apply also constraints from T{1S) — )■ Ai'j 
|61j as implemented in NMSSMTools 2.2.0. This excludes many points with a light Ai in the 
range 4 GeV < rriAi < 9 GeV. However, since this observable does not exclude a distinct 
region in (m/^±, tan/3), the constraints are not shown with a separate colour in the figures 
(though points termed 'allowed' always fulfil this constraint). We refrain from applying the 
new constraints from T(35') — ?■ Ai'j decay [U^ ESI [12] since theoretical uncertainties in the 
calculation of the exclusive Upsilon decays are not well under control [6^. 

In Figure [6] we show a comparison between the flavour constraints in scenario (A) (to the 
left) for the MSSM (upper row) and the NMSSM (lower row). We find it most useful to display 
the exclusion in the plane {mH±, tan/3). We notice the large similarities between the MSSM 
and the NMSSM allowed regions. Indeed, even if the Bg — branching ratio can be 
different in the NMSSM, the allowed region keeps a similar shape, since the extra parameters 
of the NMSSM enable many of the points to evade the Bg — )■ /^'''/^~ constraint. Allowed points 
are found down to mH± ~ 95 GeV (100 GeV) for the MSSM (NMSSM). The corresponding 
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Figure 7: Combined constraints from collider and flavour physics in the NMSSM for scenario 
(A) (left) and scenario (B) (right). The constraints are applied in the order indicated by the 
legend, with the allowed points on top. 



results for scenario (B) are given in the right column of Figure |6} The allowed region is again 
similar in both the MSSM and the NMSSM, and is larger than in scenario (A) since the lower 
limit on 5 — )■ Xg'y does not come into play. On the other hand we note that the limit on rriHi 
is stronger, mH± ^ 145 GeV in both cases. We also note that the Bg — )■ branching ratio 

is less constraining in scenario (B). 

In Figure [7| we combine the flavour constraints with the constraints from colliders presented 
in the previous section. When all constraints are applied, we do no longer observe any accepted 
points with m//± < 120 GeV. 

Knowing that the constraints from B — )■ Xsj and Bg fJ'~^f^~ depend on the choice of 
SUSY scenario (as can clearly be seen by comparing the results for our two scenarios), it is 
also motivated to look at the situation when these constraints are removed. This is presented 
in Figure [8} which only includes the direct constraints from colliders and those flavour physics 
constraints which are mediated by at tree-level. As already discussed, these are robust 
with respect to a change of the SUSY scenario. When only these 'unavoidable' constraints are 
applied, the situation at low tan/3 < 10 is once again determined by the collider constraints. 
Hence the results from the previous section applies in this region, and we flnd points with 
mH± ^ 100 GeV around tan/3 = 5. The high tan/3, low mj^i region on the other hand is 
strongly disfavoured by flavour constraints (both in the MSSM and the NMSSM). 

One constraint which was not discussed so far is that from the anomalous magnetic moment 
of the muon, = {g^ — 2)/2. As is well-known, the SM calculation of this quantity disagrees 
with the measurement of the E821 experiment. The discrepancy currently amounts to a 3.6 a 
effect [65] when hadronic contributions are evaluated with data from e~^e~ collisions, corre- 
sponding to a^^P — = (28.7 ± 8.0) x 10~^°. A slightly smaller deviation is observed when 
r decay data is used instead. A positive contribution from supersymmetry (most importantly 
from chargino/sneutrino or neutralino/smuon loops) can explain this deviation. However, this 
requires that the sign of the /i parameter is positive, which in the NMSSM translates into the 
requirement Hes > 0. Even with this choice, the SUSY mass scale in the benchmark scenarios 
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Figure 8: Combined constraints from collider and flavour physics in the NMSSM for scenario 
(A) (left) and scenario (B) (right). Only flavour observables with tree-level if^ exchange are 
included. The constraints are applied in the order indicated by the legend, with the allowed 
points on top. 

employed here is too high to explain the full deviation, but the difference between data and the- 
ory is reduced. It would be possible to explain the full difference by departing from scalar mass 
universality and allow for a lower smuon mass. This would not have an important impact on 
any of the other processes we consider. The numerical differences we obtain in between the 
MSSM and the NMSSM is typically below 1 x 10~^*^, hence much smaller than the theoretical 
error and therefore unlikely to be interesting. 

5 LHC prospects 

Given the higher centre of mass energy of the LHC compared to the Tevatron, the prospects 
for discovering NMSSM Higgs bosons is certainly much brighter at the LHC, but still remains 
a difficult task. In the following we will highlight some special features of the scenarios consid- 
ered here that are relevant for the searches at the LHC. We discuss how they differ from the 
corresponding MSSM scenarios, as well as the impact of collider and flavour constraints. We 
will concentrate on the design energy for the LHC of 14 TeV, although many of the results we 
present are equally valid for lower centre of mass energies such as 7 or 8 TeV. 

We start by considering the mass spectra of the two scenarios. As already alluded to the 
Higgs masses can differ substantially from the corresponding MSSM scenarios due to mixing 
with the additional singlet, as well as loosening of the collider and flavour constraints as dis- 
cussed in the previous section. As a main result the lightest CP-even (and CP-odd) Higgs, 
as well as the charged Higgs, can be lighter than in the corresponding MSSM scenarios. The 
relation of their masses to those of the remaining states A2, -^2, and i^s in general depends 
on which states originate from the flrst doublet, the second doublet, and from the singlet. 
The CP-odd Higgses can be classifled as either doublet-like or singlet-like, while the CP-even 
doublet states can be further separated into SM-like or MSSM-like. 
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Table 3: Reduced Higgs couplings in the NMSSM compared to the MSSM and the SM (when 
applicable). The couplings are defined such that the cross section for single Higgs and associated 
AiHj production can be obtained by multiplying the relevant cross section in Figure 10 with 
the reduced coupling squared. It should also be noted that the reduced couplings to fermions 
are identical for all three generations, even if only the third generation is displayed here. 



The CP-odd Higgs masses are related to the doublet mass m\ as < m\ and > m\, 
where the near equality holds for a purely doublet state. Through Equation ^ they are 
similarly related to w^^± • The mass of the singlet-like component is limited by Vg which varies 
inversely with A. The CP-even Higgs spectrum is slightly more complicated. Recall that in 
the MSSM there is one SM-like state with tree-level mass rrih < mz and one MSSM-like with 
mass m^jj ~ m\ = rn'jj± — m^. In the NMSSM, the mixing with the singlet state (with mass 
ms ^ Vs) will then give three states with masses m^^^ > max (m^, mj^), min (m/i, ttih) ^ ''^'h2 ~ 
max (rnh,mH), and m'jj_^ < min {mh,mH)- The (upper) equalities again holds for the pure SM 
or MSSM-like states respectively, whereas the singlet-like state can have any mass. 

In summary, the charged Higgs mass determines the overall features of the mass spectrum 
for the neutral Higgses through Equation ([s]). For large mjj± > f , both Ai and H2 will have 
masses that are smaller than or similar to the charged Higgs mass, whereas A2 and H3 will 
be heavier. Still, both A2 and can have masses of phenomenological interest for the LHC. 
However, since H3 does not add any essential new phenomena compared to H2, we will not 
discuss it further in the following. 

Using the conventions introduced in Section |2] the couplings of the Higgs bosons to fermions 
and vector bosons, as well as the triple Higgs boson couplings, can be written down. In 



19 



Table[3]we give these couplings in their reduced form, meaning that common couphng constants, 
kinematic factors etc. are left out and only the dependence on the mixing angles is kept. Thus 
the reduced couplings give a measure of the relative couplings of the NMSSM compared to 
the MSSM and the SM (where applicable). As is clear from Table |3] the NMSSM couplings 
approach their MSSM expressions in the limit cos6'a — ^ 1, cos ^13 — )■ 1, and cos 6*23 — )■ 1. We 
therefore display these mixing angles in Figure |9]for the NMSSM scans as a function of the input 
parameter mH±. The colour coding corresponds to exclusion by different constraints, which are 
applied in the order indicated by the legend. This means that the allowed points are displayed 
on top. Given the sensitivity to the SUSY scenario of BR(B — )■ Xs7) and BR(Bs — )■ 
these constraints are applied first. From the figure, we note that for small mH± the points 
in parameter space that pass the collider and flavour constraints have small mixing with the 
singlet states, meaning that the couplings are MSSM-like, whereas for larger mH± the mixing 
with the singlet states, sin 6'^ and sin6'23, can be both large and small. The same trends are 
visible in both scenarios, although this is more clear in scenario (A) since more points pass the 
constraints in this case. 



5.1 Production of neutral Higgs bosons 

The cross sections for the dominant production modes of the neutral NMSSM Higgs bosons 
can be obtained by combining the reduced couplings given in Table [3] with the production cross 
sections for SM Higgs bosons. In addition to the reduced couplings defined in Table [3] we also 
need the reduced (^gg couplings. These are defined such that the NMSSM cross sections for 
$ production can be obtained from the SM cross sections by multiplying with the reduced 
coupling squared, thus including the contributions of sparticles appearing in the loop. 



For completeness we give some of the most relevant SM cross sectionqjin Figure 10 together 
with the leading order cross section]^ for the associated AiHj production process pp ^ Z* ^ 
AiHj (the latter is given for the two cases m^. = 10 and m^. = 100 GeV). 

In addition to the modes displayed in the figure the neutral Higgs bosons can also be 
produced together with heavy quarks through processes such as 66 — )■ $ and gg — )■ tt^ with 
SM cross sections of similar magnitude as the Higgs-strahlung ones {qq — > HV). In principle 
these cross sections can be enhanced with a factor ~ tan /3 compared to the SM. However, in 
the following we will concentrate on the dominant production processes and those which are of 
special interest for a light Ai and/or and at the same time give charged leptons for signal 
tagging. 



As is clear from Figure 10 the gg production mode dominates for SM Higgs masses above 



100 GeV. In addition, since the reduced couplings HiVV are always smaller than one, the same 



^The SM cross sections in Figure 10 have been obtained from the compilation for the TeV4LHC workshop 
which can be found at http : //maltoni . web . cern . ch/maltoni/TeV4LHC/SM . h tml] In short the gg ^ H cross 
section is based on the NNLO calculation in the large top-mass limit [66l El] including soft-gluon resummation 
to NNLL [BSj- The vector boson fusion cross section has been calculated to NLO [BSJ EOl EI] using the MCFM 
program [http : // mcfm. f nal . gov and the Higgs-strahlung cross section is calculated to NNLO in QCD [75] and 
NLO in EW [75]. For more details on the calculations we refer to the website given above. 

^This cross section has been obtained using MadGraph/MadEvent [73] with the CTEQ6L1 parton densities [7S] 
and using y/s as factorisation scale. 
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Figure 9: The mixing angles cos 6a, cos and cos 623 shown as a function of the charged Higgs 
boson mass for scenario (A) (left) and scenario (B) (right). The constraints are applied in the 
order indicated by the legend. 
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Figure 10: The cross section for various Higgs production processes at the LHC {^/s = 14 TeV). 
The corresponding NMSSM cross sections can be obtained by muhiplying with the reduced 
couphngs squared. 



will also be true in the NMSSM as long as the reduced ^gg (with $ = Hi, Ai 



not too small. The reduced couplings of the relevant Higgs states are shown in Figure 11 



couplings are 
We 



leave out as mentioned above, and also Ai since it is difficult to use single Ai production 
experimentally when Ai is lightj^ Since the reduced ^gg couplings are dominated by the $66 
couplings for high tan (], the reduced couplings can in principle be of this order. However, after 
applying the flavour constraints (chiefly B — )■ rz/,- and B — )■ Dtv^), these couplings are close to 
one in both scenarios for Hi and if2; whereas for A2 they are typically smaller than one due to 
the extra suppression by sin^^- 

As has been already mentioned, the cross sections for production of neutral Higgs bosons 
through processes with electroweak bosons are typically much smaller than the ones through 
gluon fusion. Since both C{HiVV) ands C{Higg) are close to one, both production modes are 
SM-like, so the Hi — > AiAi decays may degrade the possibilities of discovery in the standard 
mode which is already very difficult in the relevant rriH^ range. In addition, using the gg — t- 
Hi — 7- AiAi channel is presumably difficult experimentally, whereas vector boson fusion is 
more hopeful [9l[77j. The electroweak production may be of phenomenological interest since it 
provides additional signatures for tagging (backward/forward jets in the case of vector bosons 
fusion and additional electroweak gauge bosons in the case of Higgs-strahlung) . 

Vector boson fusion and Higgs-strahlung production of H2 are only relevant for rriHi ^ 
150 GeV in scenario (A) since C{H2VV) is too small for large masses after applying the collider 
constraints. In scenario (B) it is almost never relevant, although this depends mainly on the 
exclusion by 6 — )• S7. This also means that the decays H2 ^ AiZ are never relevant for this 
production mode (see Figure 14). 

In the case of a light Ai, the cross section for associated AiHi production may be comparable 
to production with a vector boson, in particular for small Hi masses as shown in Figure [To] 



^There could be some hope to extract a signal in gg Aibb using the b jets to tag the event [75] . 
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Figure 11: The reduced couplings A2gg (upper), Higg (middle), and H2gg (lower) in the 
NMSSM relative to the SM for scenario (A) (left) and scenario (B) (right). The constraints are 
applied in the order indicated by the legend. 
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Figure 12: The reduced couplings AiHiZ (upper), HiVV (middle), and H2VV (lower) in the 
NMSSM relative to the SM for scenario (A) (left) and scenario (B) (right). The constraints are 
applied in the order indicated by the legend. 
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Figure 13: The relation between the HiVV and H2VV couphngs in the NMSSM. Scenario (A) 
left, scenario (B) right. 



The relevant reduced couplings for these cases are displayed in Figure 12 As can be seen in 



the figure the HiVV coupling is typically large in both scenarios after applying the flavour 
constraints. The opposite is true for the H2VV and AiHiZ couplings in scenario (A), except 
for small H2 and Ai masses respectively where again the reduced couplings can be large. 

We would like to emphasise that even though the HiVV and H2VV couplings in principle 
are smaller in the NMSSM than in the MSSM due to mixing with the singlet state, the sum rule 
'^C{HiVVy = 1 is often saturated by the two lightest states (similarly to the MSSM). This 



holds especially when flavour and collider constraints are applied, as illustrated in Figure 13 



5.2 Decays of neutral Higgs bosons 

We turn now to the decays of the neutral Higgs bosons. Again we concentrate on those modes 
which are dominant or of special interest when Ai and/or are light, still keeping in mind 
the need for charged leptons for tagging. In the case of Higgs boson production through 
gluon fusion the decay modes of most interest to us are H2 — ?■ AiZ (the branching ratio for 
Hi — )• AiZ is small, always at the level 10~^ or less) and possibly Hi — )• Ai^i, as well as 
A2 — )■ H'^W — )■ AiW^W'^ or possibly A2 — )■ HiZ — )■ AiAiZ. The same decay modes are 
also of interest for the production through electroweak gauge bosons. We note in passing that 
similarly to the MSSM, BR(Ai — t- t^t^) 0.1 above bb threshold with Ai — )■ bb saturating the 
width. Below the bb threshold BR(Ai r+r") ~ 0.9. 



Figure [14] shows that the branching ratio H2 — ?■ AiZ can be dominant. Together with the 
large cross section for H2 production this makes an interesting channel for the LHC The same 
is true for the decay A2 — ?■ H^W, whereas the branching ratio for A2 ^ HiZ turns out to be 
much smaller. This is mainly due to the comparatively smaller number of final state particles 
(not shown). Finally, we see that the Hi — )■ AiAi decay can be very dominant as soon as the 
masses are such that the decay is open. These conclusions depend only weakly on the SUSY 



scenario, as is clear from the small differences observed in Figure 14 Combining production 



and decay, the product C{A2ggY x BR(y42 — )• HiZ) is in the best case of the order 10 ^. The 
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Figure 14: The branching ratios for H2 A^Z (upper), — )■ H^W (middle), and Hi — )■ AiAi 
(lower) in the NMSSM. Scenario (A) left, scenario (B) right. 
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Figure 15: The cross sections for various charged Higgs production processes at the LHC 
{^/s = 14 TeV). The values for associated production H"^^ should be multiplied with the 
reduced couplings squared to obtain the NMSSM results. The solid line is for scenario (A) and 
the dashed one for scenario (B). 



same is true for C{A2gg)'^ x BR(A2 — ?■ H^W^), which could even so be relevant for charged 
Higgs production since this mode is of a similar magnitude as the other production modes 
(to be discussed next). 



5.3 Production and decay of charged Higgs bosons 

The main production mechanism for heavy charged Higgs bosons {mij± > rrit) at the LHC is 
through the twin processes gb tH^ and gg — > tbH^, depending on whether one is using a four 
or a five flavour scheme for the parton densities. In addition the gg — ?■ tbH^ process also gives 
the dominant production mechanism for light charged Higgs bosons {mH± < uit) which proceeds 
through top-pair production with one of the tops decaying to a charged Higgs. In both the 
NMSSM and the MSSM the magnitude of the cross section is determined by the H^tb coupling, 
gH±tb oc [nit^—^ cot (5 + mi,^-^ tan (5) , giving a minimum around tan/3 = ^Jmt/mij ~ 7. The 
production cross section above the ti threshold has been calculated to NLO both in the five 



and four ^79j flavour schemes. In Figure 15 we show the cross section calculated in the five 
flavour scheme as given by the parametrisation contained in FeynHiggs [23]. Below the ti 
threshold we have used the t -> bH^ branching ratios calculated with the same program, based 
on the NLO calculation [50], and multiplied with a ti production cross section of 900 pb jST] . 

In addition to this standard production mechanism, the possibility of a light Ai (and possibly 
a light Hi) in the NMSSM means that associated (with $ = Ai, Hi) production can also 

be large enough to be of phenomeno logical interest. This is explored in some detail in [14] 
(without the explicit application of flavour constraints). Figure 15 gives the LO cross sectioiij^ 

^Similarly to the associated HA production the cross section has been obtained using 

MadGraph/MadEvent [73] with the CTEQ6L1 parton densities [T^ and y^s as factorisation scale. For an NLO 
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m^. (GeV) m^. (GeV) 

Figure 16: The reduced HiH^W'^ coupling in the NMSSM for scenario (A) (left) and scenario 
(B) (right). 



apart from the reduced couplings HiH^W~ and AiH^W~ , which are given by Figures 16 
and |9| respectively (note that the reduced AiH^W~ coupling is equal to cosOa)- From this 
it is clear that in scenario (A) both production modes can be of interest, at least for light 
charged Higgs bosons, whereas in scenario (B) the cross section for pp — )■ HiH^ is probably 
too small to be of phenomenological interest. The dampening of C{HiH'^W~) for large mH± 
is a consequence of the sum rule C{HiVVf + C{HiH+W-f + Sf^= 1 [H] (cf. Figure |l2[. 

Finally the presence of a light Ai also opens up the possibility of a substantial branching 
ratio for — )■ AiW^ decays. As is shown in Figure 17, this remains the case in both scenarios 
also after the flavour constraints have been applied. 



6 Summary and conclusions 

To summarize, we have considered NMSSM extensions of well-known (and widely used) MSSM 
benchmark scenarios for Higgs physics. This provides a 'fair' method to compare the two 
models without resorting to theoretical bias (such as fine-tuning measures). Our numerical 
results are presented for two specific benchmark scenarios, but the method and conclusions are 
easily generalized to other scenarios. Using data from LEP, the Tevatron, and flavour physics 
experiments, we determine the allowed regions of parameter space for the Higgs bosons masses. 
In both benchmark scenarios we find that the collider limits allow for one or more Higgs bosons 
which are substantially lighter in the NMSSM {mH^ < mz, niAi ^ 10 GeV, ~ 90 GeV) 
compared to the limits that apply in the MSSM (m/, > mz, tua > mz, 'rnH± > 120 GeV). We 
also find that the low tan/3 region excluded in the MSSM is reopened in the NMSSM. These 
conclusions hold also when including present limits from channels geared specifically towards 
detection of light NMSSM Higgs bosons, such as Hi AiAi At and iy=^y4i. 

Direct limits from searches for charged Higgs bosons could in principle be more general than 

calculation and analysis of this process in the MSSM, see [52] . 
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Figure 17: The branching ratio for — )■ AiW^ in the NMSSM (computed for off-shell W^) 
in scenario (A) (left) and scenario (B) (right). 

those from neutral Higgs searches, since no new charged degrees of freedom are introduced in 
the NMSSM. Unfortunately the current sensitivity is not at the level necessary to provide any 
significant exclusion beyond the LEP results of mH± ^ 80 GeV, except for very high (and 
very low) tan /3. Improved exclusion limits on would be useful to constrain the NMSSM 
parameter space further. In fact, not even the present MSSM limits on are completely 
unavoidable, since the additional decay channel — )■ ly^Ai could degrade the performance 
of the standard searches. We therefore strongly recommend that this channel is considered for 
LHC analyses in order to fill this gap. The preliminary CDF results on this mode should serve 
as encouragement and a challenge to the LHC collaborations. 

Flavour physics experiments are complementary to those at high-energy colliders, since the 
observables measured there can be exploited to constrain the Higgs sector indirectly far beyond 
the energy scale which is directly accessible. With constraints from flavour physics included, 
we still find that lighter charged Higgs bosons than in the MSSM can be accommodated in 
the NMSSM. However, this typically requires that the underlying SUSY scenario is tuned to 
cancel the contributions to -B — ?■ transitions. Another observable of importance in the 
NMSSM is the — t- fi^fj^" decay which is sensitive to the presence of a light Ai. This leads to 
a larger variation between the new physics contributions in the MSSM and the NMSSM, since 
a light Ai can be realized independently of mH±- In the NMSSM, we find that the constraints 
from Bs — )• fJ'~^fJ'~ can be excluding down to moderate tan/3 values. It would therefore be useful 
to extend the calculation beyond the usual high tan /3 limit to study this region in more detail. 
An experimental measurement of B^ — )■ would of course also be very welcome to further 

constrain (or confirm) physics beyond the SM. For tan /3 > 10, both Bg — )■ fi^fi~ and the 
leptonic decays of pseudoscalar mesons (e.g. Bu — )■ tu) rule out the light if^ scenario. Since 
the latter are sensitive to tree-level exchange, the derived limits on (mff±,tan/3) remain 
robust in the NMSSM (as they are also with respect to changes of MSSM scenario). 
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Having applied the collider and flavour constraints, we evaluated quantities of interest for 
LHC Higgs phenomenology. Looking at the results for the reduced couplings relevant to the 
main production modes, we find that the CP-even Higgs sector in general is similar to that 
in the MSSM. This means that gluon fusion remains the dominant mode of production and 
that the coupling of the heaviest CP-even Higgs boson to vector bosons is very small. The 
main phenomenological difference here is therefore the additional decay Hi — )■ Ai^i, which 
can be large in both scenarios also after constraints are applied. In the CP-odd sector, the 
production gg — )■ A2 is suppressed due to singlet mixing. Nevertheless we find that the mode 
gg ^ A2 H^W^ deserves further investigation as a channel not present in the MSSM. 
For the direct production of charged Higgs bosons, an interesting alternative to the standard 
mode is offered by pp — >■ H^Ai, in particular when a light Ai is mostly doublet which gives a 
reduced coupling of 0{1). The same coupling also governs the decay — > W^Ai which is 
often dominant when kinematically allowed. 

In conclusion the LHC could be facing scenarios which differ substantially from the MSSM 
in the Higgs sector, while the sparticle sector is effectively unchanged. It is important to 
be aware of this possibility when designing Higgs search strategies; either trying to maintain 
model independence or by including modes which are not relevant in the MSSM. This is of 
course primarily motivated by the fact that Higgs discovery can not be guaranteed in the 
NMSSM (unlike the MSSM). Taking a more optimistic perspective, we would like to stress the 
importance that no channel is left behind since the NMSSM can offer a much richer accessible 
Higgs phenomenology that deserves to be fully exploited. 
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